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ABSTRACT 

After cosmological recombination, the primordial hydrogen gas decoupled from the 
cosmic microwave background (CMB) and fell into the gravitational potential wells of 
the dark matter. The neutral hydrogen imprinted acoustic oscillations on the pattern of 
brightness fluctuations due to its redshifted 21cm absorption of the CMB. Unlike CMB 
temperature fluctuations which probe the power spectrum at cosmic recombination, 
we show that observations of the 21cm fluctuations at z ~ 20-200 can measure four 
separate fluctuation modes (with a fifth mode requiring very high precision), thus 
providing a unique probe of the geometry and composition of the universe. 
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1 INTRODUCTION 

Resonant absorption by neutral hydrogen at its spin-flip 
21cm transition can be used to map its three-dimensional 
distribution at ea r ly cosmic times llHogan fc Reed Il979t 
IScott fc Rees!ll990l : iMadau et alJll997ft . The primordial in- 
homogeneities in the cosmic gas induced variations in the 
optical depth for absorption of the cosmic microwave back- 
ground (CMB) at the redshifted 21cm wavelength. Absorp- 
. tion occurs as long as the spin temperature of hydrogen, T„ 
(which characterizes the population ratio of the upper and 
lower states of the 21cm transition), is lower than the CMB 
temperature, T 7 . This c ondition is satisfied in the redshift 
interval 20 < z < 200 jLoeb fc Zaldarriaeal I2004T) . before 
the fi rst galaxies formed in the universe (Barkana fc Loebl 
1200 lib 

The formation of hydrogen through the recombina- 
tion of protons and electrons at z ~ 10 3 decoupled the 
cosmic gas from its mechanical drag on the CMB. Sub- 
sequently, the neutral gas was free to fall into the gravi- 
tational potential wells of the dark matter. The infalling 
gas retained some memory of the oscillations imprinted in 
the power spectrum of its density fluctuations on the scale 
traversed by sound waves in the photon-baryon fluid prior 
to decoupling. The signatures of these acoustic oscillations 
have bee n detected recently in the CMB anisotropies from 
z ~ 10 3 iSpergel et al 1l2003Tl and in t he large-scale distribu- 
tion of massive g alaxies at z ~ 0.35 fcisenstein et al]|2005t 
ICole et ail l2005f) . Here we calculate their imprint on the 



21cm brightness fluctuations at the intermediate range of 
20 < z < 200 and propose it as a probe of the composi- 
tion and geometry of the universe. Previous calculations of 
the fluctuations at these redshifts were inaccurate since they 
ignored fluctuations in the sound speed of the cosmic gas. 
Furthermore, we show that the 21cm power spectrum as a 
function of redshift can be decomposed into a combination of 
five fixed fluctuation modes, with coefficients that vary with 
redshift. We show that these modes and their coefficients 
can all be measured directly and then used to constrain cos- 
mological parameters. Several groups are currently engaged 
in the construction of low-frequency radio experiments that 
will attempt to detect the diffuse 21cm radiation from z > 6 
( http : // space .mit .edu /eor- workshop / ) . 



2 GROWTH OF DENSITY PERTURBATIONS 

On large scales, the dark matter (dm) and the baryons (b) 
are affected only by their combined gravity. The evolution of 
sub-horizon linear perturbations is described by tw o coupled 
second-order differential equations ([Peebles 1983): 



5dm + 2H6 dm 

S h + 2H5 b 



47rGp m (/b5b + /dm^dn 
4nGp m (fb$b + /dm<$dn 



(1) 
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where Sdm(t) and Sb(t) are the perturbations in the dark 
matter and baryons, respectively, the derivatives are with 
respect to cosmic time t, H(t) = a/a is the Hubble constant 
with a — (1 + z)" 1 , and we assume that the mean mass 
density p m (t) is made up of respective mass fractions /dm 
and /b = 1 — /dm- Since these linear equations contain no 
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spatial gradients, they can be solved spatially for 5d m (x,i) 
and <5b(x,i) or in Fourier space for 5d m (k,t) and 5b(k, t). 
Denning 5tot = fbfib + /dm^dm and 5b- = 8b — 5tot , we 

find 



<5tot + 2HStot 
5 b - + 2if5 b _ 



. 



(2) 



Each of these equations has two independent solutions. The 
equation for 5tot has the usual growing and decaying so- 
lutions, which we denote D\(t) and D4,(t), respectively, 
while the 5b- equation has solutions D^it) and Ds(t); we 
number the solutions in order of declining growth rate (or 
increasing decay rate). We assume an Einstein-de Sitter, 
matter-dominated universe in the redshift range z = 20- 
150, since the radiation contributes less than a few percent 
at 2 < 150, while the cosmological constant and the curva- 
ture contribute to the energy density less than a few per- 
cent at z > 3. In this regime a oc t 2 ^ 3 and the solutions are 
D±(t) = a/ai and D 4 (t) = (a/a;)~ 3/2 for Stat, and D 2 (t) = 1 
and D3(t) = (a/ai)^ 1 ^ 2 for 5b_, where we have normalized 
each solution to unity at the starting scale factor a^. Intu- 
itively, since baryons and dark matter both feel the same 
gravity, the difference between them changes according to 
D2 and D3 which describe the evolution of perturbations in 
the absence of gravity. We set the initial redshift Zi = 150, 
when the photon density and temperature perturbations be- 
come negligible o n sub- horizon scales and can be neglected 
l|Naoz fc Barkanall2005l) . 

The observable baryon perturbation can then be written 

as 

4 

5 b (k,t) = 5 b - + Stot = M k ) D ™W > ( 3 ) 

m— 1 

and similarly for the dark matter perturbation, 

4 

5dm = (5 tot - f b 5 b ) = V <L(k) C m (t) , (4) 

/dm 

m — 1 

where d = Di for i = 1,4 and d = -(/b//dm)A for i = 
2, 3. We establish the values of 5 m (k) by inverting the 4x4 
matrix A that relates the 4- vector (5i, 82, 53, 54) to the 4- 

vector that represents the initial conditions (5b, 5d m , 5b, 5d m ) 
at the initial time. 



where or is the Thomson cross-section, x e (t) is the elec- 
tron fraction out of the total number density of gas par- 
ticles, and p 1 is the CMB energy density at a tempera- 
ture T 7 . In the redshift range of interest, we assume that 
the photon temperature (T 7 = T"/a) is spatially uniform, 
since the high sound speed of the photons (i.e., c/s/3) sup- 
presses fluctuations on the sub-horizon scales that we con- 
sider, and the horizon-scale ~ 10 -5 fluctuations imprinted 
at cosmic recombination are also negligible compared to the 
smaller-scale fluctuations in the gas density and tempera- 
ture. Fluctuations in the residual electron fraction x e (t) are 
even smaller. Thus, 



^ _2JlWb 4), „, -4 

~dt~3 1 ^^ + ~tT {l1 I)a 



(7) 



where t' 1 = p° (8cr T c/3m e ) = 8.55 x 10" 13 yr _1 . After cos- 
mic recombination, x e (t) changes due to the slow recombi- 
nation rate of the residual ions: 



dx e (t) 
dt 



= -a B (T)x 2 (t)n H (l + y) , 



(8) 



where as(T) is the case-B recombination coefficient of hy- 
drogen, TiH is the mean number density of hydrogen at time 
t, and y = 0.079 is the helium to hydrogen number density 
ratio. This yields the evolution of the mean temperature, 
df/ dt = -2HT + x e (t)t~ 1 (T-y -T) a~ 4 . In prior analyses 
(e.g.. lPeeblesil980l : lMa fc Bertschingerll995D a spatially uni- 
form speed of sound was assumed for the gas at each red- 
shift. Note that we refer to Sp/8p as the square of the sound 
speed of the fluid, where Sp is the pressure perturbation, 
although we are analyzing perturbations driven by gravity 
rather than sound waves driven by pressure gradients. 

Instead of assuming a uniform sound speed, we find the 
first-order perturbation equation, 



dSr 
~df 



2 dS b 
3~dT 



(9) 



where we defined the fractional temperature perturbation 
St- Like the density perturbation equations, this equation 
can be solved separately at each x or at each k. Furthermore, 
the solution 5y(i) is a linear functional of Sb(t) [for a fixed 
function x e (t)]. Thus, if we choose an initial time U then 
using Eq. © we can write the solution in Fourier space as 



5 T (M) = ^5 m (k)L>^(t)+5 T (k,t i )i?J(i) 



(10) 



3 GROWTH OF TEMPERATURE 
PERTURBATIONS 

We include the (previously ignored) fluctuations in the 
sound speed of the cosmic gas caused by Compton heating 
of the gas, which is due to scattering of the residual electrons 
with the CMB photons. The evolution of the temperature 
T of a gas element of density pb is given by the first law of 
thermodynamics: 



dQ = 2 fc srfT - fc B Tdlogp b 



(•») 



where dQ is the heating rate per particle. Before the first 
galaxies formed, 



^ =4^1k B {T 1 -T)p^x e {t) 
at m e 



(6) 



where D^(t) is the solution of Eq. © with St = at 
ti and with the perturbation mode D m (t) substituted for 
5f,(t), while Dq{£) is the solution with no perturbation S b (t) 
and with 5 T = 1 at U By modifying the CMBFAST code 
iSeliak fc Zaldarriagalll99al . we numerically solve Eq. JUJ 
along with the density perturbation equations for each k 
down to Zi = 150, and th en match the solution to the form 
of Eq. UHI . We note that lBharadwai fc A\i (l2004al) derived 
a similar equation to Eq. © but solved it only for a den- 
sity perturbation that follows the growing mode D\(t), ne- 
glecting the other density modes and the initial temperature 
perturbation. 

Figure shows the time evolution of the various in- 
dependent modes that make up the perturbations of den- 
sity and temperature, starting at the time ti correspond- 
ing to Zi = 150. D^{t) is identically zero since D2(t) — 1 
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Figure 1. Redshift evolution of the amplitudes of the indepen- 
dent modes of the density perturbations (upper panel) and the 
temperature perturbations (lower panel), starting at redshift 150. 
We show m = 1 (solid curves), m = 2 (short-dashed curves), 
m = 3 (long-dashed curves), m = 4 (dotted curves), and m = 
(dot-dashed curve). Note that the lower panel shows one plus the 
mode amplitude. 



is constant, while Dj(t) and Dj(t) are negative. Figure[5] 
shows the amplitudes of the various components of the ini- 
tial perturbations. We consider comoving wavevectors k in 
the range 0.01 - 40 Mpc -1 , where the lower limit is set by 
considering sub-horizon scales at z = 150 for which photon 
perturbations are negligible compared to 8dm and <5b, and the 
upper limit is set by requiring baryonic pressure to be neg- 
ligible compared to gravity. 82 and £3 clearly show a strong 
signature of the large-scale baryonic oscillations, left over 
from the era of the photon-baryon fluid before recombina- 
tion, while Si, 84, and St carry only a weak sign of the oscil- 
lations. For each quantity, the plot shows [fc 3 P(k) / (2n 2 )] 1 ' 2 , 
where P(k) is the corresponding power spectrum of fluctu- 
ations. 

84 is already a very small correction at z = 150 and 
declines quickly at lower redshift, but the other three modes 
all contribute significantly to Sh, and the St^U) term remains 
significant in <5t(£) even at z < 100. Note that at z = 150 
the temperature perturbation St has a different shape with 
respect to k than the baryon perturbation 8b, showing that 
their ratio cannot be described by a scale-independent speed 
of sound. 

The power spectra of the various perturbation modes 
and of <5r(ii) depend on the initial power spectrum of den- 
sity fluctuations from inflation and on the values of the fun- 
damental cosmological parameters (S7dm, ^6, ^a, and h). If 
these independent power spectra can be measured through 
21cm fluctuations, this will probe the basic cosmological 
parameters through multiple combinations, allowing consis- 
tency checks that can be used to verify the adiabatic nature 
and the expected history of the perturbations. Figure 



Figure 2. Power spectra and initial perturbation amplitudes ver- 
sus wavenumber. The upper panel shows <5t, (solid curves) and 
5d m (dashed curves) at z = 150 and 20 (from bottom to top). 
The lower panel shows the initial (2 = 150) amplitudes of Si 
(solid curve), 62 (short-dashed curve), £3 (long-dashed curve), 84 
(dotted curve), and 5t(*j) (dot-dashed curve). Note that if 5i is 
positive then so are £3 and <5r(ti), while §2 is negative at all k, 
and 64 is negative at the lowest k but is positive at k > 0.017 
Mpc -1 . 



lustrates the relative sensitivity of P(k) to variations in 
Sljm/i 2 , Qbh 2 , and h, for the quantities Si, 82, 83, and <5t(£i)- 
Not shown is 84, which although it is more sensitive (chang- 
ing by order unity due to 10% variations in the parameters) , 
its magnitude always remains much smaller than the other 
modes, making it much harder to detect. Note that although 
the angular scale of the baryon oscillations constrains also 
the history of dark energy through the angular diameter 
distance, we have focused here on other cosmological pa- 
rameters, since the contribution of dark energy relative to 
matter becomes negligible at high redshift. 



4 21CM FLUCTUATIONS 

The spin temperature T s is defined through the ratio be- 
tween the number densities of hydrogen atoms, ni/no = 
(31/50) ex P {—T*/T s } , where subscripts 1 and correspond 
to the excited and ground state levels of the 21cm transition, 
(31/50) — 3 is the ratio of the spin degeneracy factors of the 
levels, and T* = 0.0682K corresponds to the energy differ- 
ence between the levels. The 21cm spin temperature is on the 
one hand radiatively coupled to the CMB temperature, and 
on the other hand couple d to the kinetic gas temperature T 
through collisions (Allison & Dalgarno 1969 ,) or th e absorp- 
tion of Lya photons Twou^iuvBenl ll952l:lFieldll"95i). For the 
conco rdance set of cosmological parameters (ISpergel et all 
2003), the mean brightness temperatur e on the sky at red - 
shift z (relative to the CMB itself) is jMadau et al.lll997l) 
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Figure 3. Relative sensitivity of perturbation amplitudes at z = 
150 to cosmological parameters. For variations in a parameter x, 
we show dlog-^/ P (k) / dlog(x) . We consider variations in fl^h 2 
(upper panel), in Slj,h 2 (middle panel), and in the Hubble constant 
h (lower panel). When we vary each parameter we fix the other 
two, and the variations are all carried out in a flat f2 tota i = 1 
universe. We show the sensitivity of Si (solid curves), &2 (short- 
dashed curves), £3 (long-dashed curves), and 5y(tj) (dot-dashed 
curves). 

T b = 28 mK [(1 + z)/10] 1/2 [{T s - T 7 )/T s ] xm, where xm is 
the mean neutral fraction of hydrogen. 

In general, fluctuations in Tb can be sourced by fluctu- 
ations in gas density, temperature, neutral fraction, radial 
velocity gradient, and Lyq flux from galaxies. The veloc- 
ity gradient term <ISobolevlll960l) is in Fourier space jKaiserl 

Il987t iBharadwai fe Alii l2004allb|) 5 drVr = -fi 2 8/H, where 
H = cos Ok in terms of the angle 6k of k with respect to 
the line of sight. We ca n therefore write the fluctuation as 
jBarkana fc Loebll2005l) 

*r b (M) =^l{\i)H~ 1 +p8(k)+p T 5 T {ii)+l^) , (11) 

where we have defined time-dependent coefficients /3 and 
Pt (combining the rel evant explicit terms from Eq. 2 of 
iBarkana fe Loebl <l2005l) ). and have also combined the terms 
that depend on the radiation fields of Lya photons and 
ionizing photons into 5 ra d- Before the first galaxies formed, 
xhi = 1 and <5 ra d = 0. Thus, 



MM) = 



4 

E 



Mk) [n 2 H-\t)D m {t)+f3D m {t) 



+/3rD^(t)] + f3 T S T (k, U) Dl(t) . (12) 



The power spectrum Pr b is defined by 
(fe i (k 1 ,t)fe ilt (k a )> = (27vfS D (k 1 +k 2 )P Tb (k 1 ,t) 



(13) 



At z — 150, the overall contributions of the Si, 62, S3, and 
5t(U) terms to Pr h are comparable, while the 84, term (which 
corresponds to the usual decaying mode) represents only a 



1% correction to the 21cm power spectrum. The contribu- 
tions of the different terms to Sr b then scale with redshift 
roughly like the corresponding density modes D m {t), with 
the Sriti) term varying roughly as oc 1/a, so that the usual 
growing mode Si dominates at low redshift. 

Noting that observations of the 21cm power spectrum 
can be analyzed as a function of k and /x, we write the p ower 
spectrum as a polynomial in iBarkana fc Loebll2005 



P Tb (k) = M 4 P M 4 (fc) + /i 2 P M 2 (fc) + P M o (fc) 



(14) 



and obtain the three separate power spectra of 21cm fluctu- 
ations. Figure^Jshows their values at various redshifts. The 
fluctuations rise with time until z ~ 50 due to the growth 
of fluctuations and then decline at lower redshift due to the 
reduced collisional cou pling of the 21cm spin tem perature 
to the gas temperature llLoeb fc ZaldarriagalEooil . The five 
power spectra shown in the lower panel of Fig. [5] along with 
the coefficients of the five modes in Eq. 11211 can all be de- 
termined directly from observations. Suppose, e.g., that we 
measure the three observable power spectra at N fc-values 
at M different redshifts, for a total of 3NM data points. As 
we have shown, these data points can be modeled in terms 
of the amplitudes of five separate fluctuation modes at the 
same N fc-values, along with two coefficients per redshift for 
each fluctuation mode [corresponding to the contributions 
to the /j 2 term and to the (1- independent term in Eq. 1121 1. 
The total of 5iV + 10M parameters can then be determined 
as long as the number of data points is larger than this total. 
Thus, M = 4 redshifts suffice as long as N ^ 5, while M = 2 
suffices for any N ^ 20. Measurements at additional red- 
shifts will allow a more accurate reconstruction along with 
multiple consistency checks. Note that Figure shows that 
separating out the different modes at z ~ 20 requires an 
order of magnitude higher measurement precision than at 
z ~ 100. 

Detection of these fluctuations with future experiments 
JZaldarriaea et alJ Eool iMorales fc Hewitti Eool iMorakd 



2005) would set new constraints on the evolution of the 
universe through a regime of cosmic time that has never 
been probed before. Unlike existing cosmological probes, the 
three-dimensional nature of 21cm fluctuations makes possi- 
ble measurements at a wide ra nge of redshifts and scales. 
CMB temperature fluctuations jspereel et alj|2005l) probe 
the power spectrum only at cosmic recombination, and only 
on scales fc ^ 0.2 Mpc - which are not damped by pho- 
ton diffusion. Ga laxy redshift surveys dTegmark et alJl2004l : 
ICole et al.1 120051) probe the density power spectrum only 
at low redshifts where the fastest-growing mode Di(t) is 
already strongly dominant, and the signature of the bary- 
onic oscillations is weak. Furthermore, even on large scales 
(fc < 0.1 Mpc -1 ) the biased distribution of galaxies may be 
affected by complex a strophysical feedback p rocesses such 
as cosmic reionization iBarkana fc Loeb 2004). 

Indeed, CMB fluctuations and galaxy surveys cannot in- 
dividually constrain multiple-parameter cosmological mod- 
els due to degeneracies, and must be combined in order to 
achieve significant constraints (iTeemark et all2004l) . In con- 
trast, 21cm fluctuations provide a signal at a sufficiently 
early epoch when the baryons have not yet caught up with 
the dark matter, and four separate fluctuation modes are de- 
tectable (with a fifth mode requiring higher precision). Since 
each mode depends differently on the cosmological param- 
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Figure 4. Power spectra of 21cm brightness fluctuations versus 
wavenumber. We show the three power spectra that are sepa- 
rately observable, P^i (upper panel), P^2 (middle panel), and 
P ' o (lower panel). In each case we show redshifts 150, 100, 50 
(solid curves, from bottom to top), 35, 25, and 20 (dashed curves, 
from top to bottom). 



eters (see Fig. such measurements would provide great 
discriminatory power in determining the cosmological pa- 
rameters and in testing for any deviations from the stan- 
dard ACDM paradigm. This cosmological information is in 
principle available from measurements in any redshift range, 
but the amplitudes of some of the modes decline at low 
redshift and they become more challenging to measure (see 
Fig-Q. While we have focused on deriving constraints from 
the power spectra, we note that the coefficients in Eq. 112H 
depend on known cosmological and atomic physics, so that 
measuring them from the 21cm data will help to further 
constrain the fundamental cosmological parameters. 
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